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Supersymmetric spectra with a stop NLSP and a neutralino or gravitino LSP present a 
special challenge for collider searches. For stop pairs directly produced from QCD, the 
visible final-state particles are identical to those of top quark pair production, giving very 
similar kinematics but often with much smaller rates. The situation is exacerbated for 
compressed spectra with m~ t ~ m t + ?tilsp, as well as for lighter stops which can suffer from 
low acceptance efficiencies. In this note, we explore the power of a direct stop search using 
dileptonic m^ 2 , similar to the one recently performed by ATLAS, but more optimized to 
cover these difficult regions of the (m^m^sp) plane. Our study accounts for the effects 
of stop chirality and LSP identity, which can be significant. In particular, our estimates 
suggest that m 4 - ~ m t with a massless LSP is excludable for right-handed stops with bino- 
like (gravitino) LSP with 2012 (2011) data, but remains largely unobservable in the case of 
a higgsino-like singlino LSP. For each of these cases we map out the regions of parameter 
space that can be excluded with 2012 data, as well as currently allowed regions that would 
yield discovery-level significance. We also comment on the prospects of a precision m,T2 
shape measurement, and consider the potential of ATLAS's dileptonic t — > bxt searches 
when re-interpreted for light stops decaying directly to the LSP. 



I. INTRODUCTION 



Weak-scale supersymmetry (SUSY) is perhaps the most well-studied extension of the 
Standard Model (SM), since it offers a weakly-coupled solution the to Higgs hierarchy 
problem, as well as a dark matter candidate and an indication of grand unification at 
high energies. However, LHC searches have rapidly been excluding "generic" superparticle 
spectra with masses below roughly 1 TeV jll, 0]. This places significant tension on super- 
symmetric models, which require in particular that stops should be lighter than about a 
TeV to avoid the re-introduction of fine-tuning. In order to preserve the main virtue of 
supersymmetry while avoiding these constraints, we are led to seriously consider the pos- 
sibility that third generation squarks are the only colored superparticles in the sub- TeV 
range [3l-lo|. 

Standard SUSY searches rely heavily on the assumption of copious production of colored 
superparticles and the release of a substantial amount of missing transverse energy ($ T ) 
in their cascade decays, ending in a stable and invisible lightest supersymmetric particle 
(LSP). In scenarios where the LHC can only access squarks of the third generation, these 
searches become significantly less sensitive. To a large extent this is due to the reduction 
in the number of colored SUSY production channels, and the fact that third generation 
squarks cannot be efficiently produced through gluino exchange. While the cross sections 
become much higher if these squarks are somewhat light, the visible and invisible energy 
released then tends to be much smaller, and the signal can easily be missed in searches 
tailored to more generic spectra. This has motivated a new interest in dedicated searches 
for light third generation squarks. 

One of the simplest and best-motivated channels for such a search is QCD pair produc- 
tion of the lightest stop eigenstate, with each stop decaying into a top quark and the LSP. 
The particles in the final state are then identical to those in top quark pair production, 
supplemented by the two invisible LSP's. If there is a large mass gap between the stop and 
the LSP, stop production events can be distinguished from SM ti production by identify- 
ing kinematically extreme events: events with large $ T in the all-hadronic channel, events 
with large transverse mass mr (l,#rr) in Z+jets channel, and events with large "stransverse 



mass" mr2^ + ,i ,$t) Hr4l3j] in the dileptonic channel. Indeed, several LHC searches are 



already capitalizing on these strategies, and are beginning to set meaningful limits 



14- 



20l | . Nonetheless, there remain regions of parameter space that are not well-covered by 



these searches, either because the missing energy release is too small (the compressed case 
mi ~ m t + ttilsp) or because the visible energy release is too small (the off-shell case 
mi < m t + ttilsp)- For example, in the case of mi ~ m t and a nearly- massless LSP, the 
final-state kinematics become very similar to SM ti production, but with cross sections that 
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are roughly six times smaller and which rapidly decrease as the stop and LSP masses are 
raised together. Lighter stops can have much larger cross sections, but also much smaller 
acceptances, and the latter can outweigh the former depending on exactly how a search is 
defined. 

In this paper, we will study the kinematics of events arising from compressed spectra 
and spectra with off-shell decays, and survey regions of the (m^, 771lsp) plane where signal 
sensitivity could be enhanced with a more optimized search strategy using dileptonic m^- 
Optimization of light stop searches is a task that has recently been taken up by several 
other groups, which used a diverse set of approaches. In [2l|, a data-driven, shape-based 
analysis of the $ T and ray distributions of all-hadronic and /+jets channels was considered. 
Good sensitivity in the compressed limit was claimed for 2012, with exclusion up to rrif ~ 
325 GeV. However, this study effectively considered statistical errors only, and not the 
possible systematic uncertainties involved in extrapolating from control-region fits to signal- 
region predictions. Nor is it clear how to define control regions that are enriched in tops 
but depleted in stops at high-/^ T or high-m T , as would be necessary to obtain a precise 
background-only shape prediction from the data. Ref. 22( also considered a high-precision 
kinematic discrimination, but instead of exploiting the different decay kinematics of stops 
versus tops, differences in production kinematics and spin correlations were used. The 
former was revealed using the rapidity distributions of the top/stop decay products, and the 
latter was revealed using a matrix-element style discriminator. In the dileptonic channel, 
these methods were shown to yield 2-3cr discrimination between SM top production and 
tops plus stops by the end of 2012, for the specific case of mi ~ 200 GeV and m^o ~ 0. 
These results were also obtained in the absence of experimental systematics, but they are 
expected to be quite stable against theory systematics. The case of stops decaying to a 



gravitino LSP in gauge mediation was considered in |23H25|. Using the results of various 
Tevatron and LHC searches (up to 1 ftr 1 at 7 TeV) not dedicated to this particular process, 
exclusion potential for mi < m t was claimed. These studies also discussed the potential 
of dedicated searches using observables such as m(l,b). Additional studies that consider 
strategies to improve s top searches, but do not specifically address compressed or off-shell 



cases, can be found in 26|-|28|. 



To this long list of preceding work, we seek to add the following. First, we will design 
a simple "cut and count" style search with the dileptonic 777,7-2 variable, optimized for 
compressed spectra, that allows for the presence of nonvanishing systematic errors. For the 
recent ATLAS searches these errors tend to be at the O(10%) level, and we will assume 
systematics of this size. As we will show, our strategy also automatically yields sensitivity 
to off-shell decays. Second, we will investigate the impact of the model-dependence that is 
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present even in this highly simplified SUSY spectrum. ATLAS searches for stops decaying 
directly to the LSP always assume a right-handed chiral stop and a bino-like neutralino. 
Most theory papers make the same assumption, or in some cases assume a gravitino LSP. 
For our search, we will also consider what happens if left-handed stops contribute, or if the 
stop decays into a light Higgsino-like neutralino without a nearby chargino, as can easily 
arise if a light singlino mixes into the Higgs sector. We will show that we can achieve very 
good exclusion coverage in some cases, while others are more difficult. We will also discuss 



how these model-dependent effects become relevant for searches in the I- 



While ATLAS has already performed a search using dileptonic m^2 Is], that search was 



■jets channel. 



optimized for finding stops well above the compressed regime, where a large amount of lfi T 
is released. Thus an my 2 cut significantly beyond my/ could still have reasonable signal 
acceptance. Currently, the lowest excluded m t ~ with vanishing neutralino mass from this 
search is 280 GeV. We propose simple modifications which greatly enhance the potential to 
exclude or discover stops at lower masses, even those in the compressed limit and below. We 
will argue that applying these changes to the ATLAS 2011 dileptonic analysis can extend 
their exclusion boundary at m^o = from 280 GeV down to 210 GeV, and simultaneously 
exclude most of the region mj < m t that is still allowed by LEP. This still leaves a gap 
in exclusion- level sensitivity between 175 and 210 GeV, but we will show that this gap 
can be completely closed using 2012 data. The gap inevitably re-opens when we consider 
more massive neutralinos. We estimate that its width is typically 5~10 GeV, and will 
shrink further when the 14 TeV LHC comes online. Our study also reveals regions in the 
(m^, iTiyo) plane with discovery-level potential near and below the compressed limit, which 
are currently not excluded. As a supplemental analysis, we re-interpret ATLAS 's light 
t — > b\i searches, which can cover direct t — > blv\\ in complementary regions where the 
stop/LSP mass splitting is too small to produce an on-shell W . 

These results apply to the nominal ATLAS model with a right-handed stop decaying 
into a bino-like neutralino. When we consider more general models, the distributions 
can change significantly due to a combination of momentum-scaling and spin effects in the 
decay matrix elements. We find that the most promising case is actually a right-handed 
stop decaying to gravitino, tn — > tG, for which we achieve excellent coverage for essentially 
arbitrary stop masses below 380 GeV using only 2011 data. If we instead look at ti — > tG, 
then the rriT2 analysis becomes less effective and a gap remains between = 190 GeV and 
240 GeV even with the full 2012 dataset. The most difficult cases are a right-handed stop 
decaying to a higgsino-like LSP — > th°), or a left-handed stop decaying to a bino-like 
LSP (ti —> tB). These are phenomenologically equivalent, and have the softest spectra 
amongst our models, in some cases even softer than ti itself. For these models, stops of 
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less than 240 GeV with a massless LSP cannot be excluded with the 2012 data using our 
method. Compressed spectra must then be covered using other precision measurements, 
such as those proposed by 



21 



22]. 



The good results obtainable with a simple cut-and-count style search suggest that dilep- 
tonic rriT2 would itself be an excellent candidate for a precision shape-based measurement. 
Neglecting systematics and again focusing on a right-handed stop decaying to bino-like 
neutralino, we estimate that stops of any mass near the compressed limit could be ruled 
out for m^o < 40 GeV by the end of 2012, and up to nearly 200 GeV after the energy 
upgrade. 

The outline of the paper is as follows. In section [III we review the kinematic variable 
7BT2 and motivate its use in our search. In section IHH we discuss the physics of stop decay 
in more detail. We show how nontrivial distortions of m^2 distributions persist even in 
the compressed limit, and explain their mo del- dependence. In section IIV[ we begin by 
reproducing the results of the ATLAS dileptonic m^2 analysis in order to validate our 
methods, and proceed to show how a modified set of cuts can greatly expand the reach for 
light stops. In section [V] we evaluate the prospects of our search strategy for all choices of 
stop chirality and LSP identity, and consider some of its limitations and possible extensions. 
We summarize our findings and conclude in section IVII Two appendices describe in more 
detail our simulations and statistical procedures. 



II. DILEPTONIC m T2 IN STOP SEARCHES 

For each of the three main ti decay channels, we can construct a simple kinematic 
observable that, under ideal circumstances, exhibits an edge or endpoint. In the case of 
all-hadronic decay, this observable is simply and the "edge" is at zero. For Z+jets decay, 
the observable is the transverse mass of the lepton-^ T system, mr(i,^). This is the mass 
obtained by projecting the entire event into the transverse plane, pretending that the lfi T 
vector arises from a single massless particle, and taking the lepton-^ T system mass in 2+1 
dimensions. In events with a single leptonic W, rriT is automatically smaller than the full 
W mass, and has a endpoint at mw This endpoint is sculpted into a sharply-peaked edge 
by a phase space Jacobian factor. 

For the dileptonic decay, the analog of my is the "stransverse mass" my 2 (I + , l~,$ T ) jllT - 



13| . built out of the subsystem consisting of the two leptons and ~$ T . To motivate the 
construction of m^, we assume that the measurable (l + ,l~ ,$ T ) system arises from two 
identical one-step decay chains X + — > l + v and X~ — > l~v. If it were somehow possible to 
find the correct decomposition of $ T into into its individual neutrino contributions, and 
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break the combinatoric ambiguity over which neutrino to assign to which lepton, we could 
estimate the mass of the hypothetical parent X by again constructing m^'s for both sides 
of the event. Of course the correct $ T decomposition is not available in practice, however 
there is a procedure to obtain a lower bound for the individual rrty's on an event-by-event 
basis. This is achieved by considering all possible partitions Pt{v) + prip) = $t-> an d 
computing the following extremum: 

rriT2{l Jr \l~ ,$t) = mrn [max{mr(/ + , v), mr(l~, 9)}] . (1) 

Pt(v)+Pt{v) = $ T 

In events which genuinely contain two on-shell W's serving as our X's, and no additional 
sources of $ T , the correct transverse v and v configuration yields mx{l + , v) and mr{l~ , v) 
that are both smaller than mw- Since the above procedure by definition gives a number 
that is less than or equal to the larger of these two true %'s, m^2 is also strictly bounded 
above by mw- The endpoint at mw is not technically an edge, since the distribution falls 
continuously to zero over a finite range, but the approach is nonetheless quite steep. 

When all final state particles are massless, there is an analytic procedure to compute 
29|. The trajectories of the leptons in the transverse plane can be viewed as defining 
the boundary of two wedges, one of which has an opening angle less than it. If the lfi T 
vector lies within this wedge, then the configuration that extremizes equation ([1]) is the one 
where the two neutrinos move collinearly to their sister leptons, and in this configuration 
mx2 is identically zero. Otherwise, the extremization problem can be reduced to finding 
the roots of a quartic equation. The extremized configuration has A<f)(v, v) = — A0(/ + , l~) 
and my(/ + ,z/) = mr(!",i/). In the special case where Jfi T exactly balances the pr of the 
two leptons, such as we would get with two W's with equal and opposite p^s (including 
two W's at rest), the solution is just (pt(v),pt(v)) = (—Pt(1~),—Pt(1 + ))- Thus in this 
special kinematic configuration: 



m T2 ->■ ^2p T {l+)p T {l-) (1 + cosA0(/+,/-)). (2) 

The largest m-^ values, near mw, are obtained when both W decays are fully transverse 
to the beampipe, and the two leptons are collinear and balanced against the two neutrinos. 
While derived for a special case, this still gives an approximate picture of the kinematics 
as long as the W + W~ system is not itself highly boosted in the transverse plane. 

A stop search can utilize any of the three final states and the corresponding observables 
mentioned above: $ T for all-hadronic, m<r for Z+jets, and mx2 for dileptonic. The two LSP's 
produced in the stop decays inject additional $ T into the event, allowing values beyond 
the edges/endpoints in each of these observables to be populated, in principle leading to 
significant enhancements in signal over background. 
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Of course, taking the realities of a hadron collider environment into consideration, the 
purification of S/B in these three channels is not equally practical. All-hadronic tt events 
have Ifirp from sources such as jet energy mismeasurement and 5-meson decays, and due 
to contamination from tt decaying into Z+jets or even dileptonic final states (especially 
those with taus). It is therefore quite common to observe all-hadronic events with 10's 
of GeV of 1fi T , which cause a significant background for genuine signal events, especially 
in compressed _spectra. Nonetheless, all-hadronic searches are possibly the best option for 
heavier stops 14lll7H20|. and it will be important to fully understand how far these searches 



might be pushed into the compressed regime. 

The next-best strategy is to conduct a search in the Z+jets channel, looking for events 
with tut > rn^. However, a search in this channel also suffers from tt backgrounds be- 
yond the kinematic edge, mainly due to misidentified dileptonic events with a hadronic 
tau. The situation is nonetheless much better than for all-hadronic, and ATLAS exclusion 
limits based on the semileptonic final state currently extend down to m- t ~ 225 GeV for a 
massless bino-like LSP and right-handed stop 15j. But this result is already systematics 
limited, and a significant improvement will likely require better modeling of the dileptonic 
contamination. Another possibility, already utilized in the all-hadronic search, is to place 
an explicit tau veto, perhaps supplemented with a soft lepton veto. 1 We do not explore 
that option here, as it is difficult for us to reliably model. 

If we are simply interested in a basic cut-and-count analysis with reasonable S/B, that 
leaves us with the dileptonic channel and the rriT2 variable. Of the three different variables 
and decay channels, this one has the most well-defined endpoint with the fewest outliers, as 
has clearly been observed experimentally [l6|. Events with one or more taus are harmless, 
as only leptonic taus can contribute, and these tend to populate lower m^ 2 values. The 
few background events with far beyond mw tend to arise either from tt events with 
mismeasured lfi T , or from SM electroweak backgrounds and tiW/Z. 

In Fig. [TJ we show fully reconstructed dileptonic m T2 and semileptonic m T distributions 
for two very moderately compressed example models (just below the current ATLAS Z+jets 
exclusion) and for the dominant background tt. The distribution uses cuts that we 
describe in section HVBt while the mx distribution uses cuts defined by ATLAS for their 
lowest-mass search region [I5J]. 2 The former clearly has a significantly better S/B. The 



1 We estimate that approximately 2/3 of ATLAS 's high-m^ events are dileptonic tt with a tau. Most of 
the rest are prompt dileptonic events with a lost lepton. 

2 We have verified that our own Z+jets analysis faithfully reproduces that of ATLAS. In particular, we 
obtain the same shape of rriT with preselection cuts, and accurately predict the tt counts in all five of 
their signal regions. 
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m T2 (GeV) m T (GeV) 

FIG. 1: The distribution of dileptonic m^2 (left) and semileptonic (right) for 220 GeV right- 
handed stops decaying into either a massless bino-like neutralino or a massless higgsino-like neu- 
tralino. For comparison, the distribution for the SM ti background is also shown. Events have 
been processed through the reconstruction described in appendix The txiti distribution was 
formed after application of our final analysis cuts described in section IIVB1 The my distribution 
was formed after application of the "SR. A" type cuts described in [3] , except for the cut on rriT 
itself. (Error bars are Monte Carlo statistics.) 

difference between the models is due to spin effects, which we will shortly discuss in detail. 

III. STOP DECAY 
A. Basic kinematics 

In a simplified SUSY spectrum containing only the lightest stop eigenstate and an in- 
visible LSP, stops decay into the same final states as top quarks, but with the addition of 
the LSP. In high-scale mediation models with a neutralino LSP, a semileptonic stop decay 
looks like t — > b^vxi- In low-scale mediation models with a gravitino LSP, the analog is 
t — > bl + uG. If the stop is heavy enough relative to the LSP, namely m t ~ > m t + mLSP, it 
can simply undergo a 2-body decay t — > t(xi/G) followed by the decay of the on-shell top 
quark. Indeed, to date, two body decays are always assumed in LHC searches for stop pairs. 
However, if the stop is somewhat lighter, but still heavy enough to produce an on-shell W 
boson {m w + m b + m LS p < m t - < m t + m LSP ), it will instead undergo a 3-body decay 
t — > bW + (xi/G). Even lighter stops, with mi < mw + m^ + ?tilsp, will naively undergo a 
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fully 4-body decay. 3 Ultimately, as the available phase space is completely closed off, the 
stop becomes quasi-stable, a case which we do not consider. 



As has been pointed out before 



24|, the transitions between the different iV-body 



decay regions are not perfectly abrupt. For example, consider a maximally compressed 
spectrum with a massless LSP and a stop infmitesimally heavier than the top quark. The 
phase space for the 2-body decay is nearly vanishing, naively providing us with an LSP 
with zero 4-momentum and a final state that looks exactly like a top quark. However, in 
reality the 3-body decay (with m(W, b) < m t ) will have much more phase space available in 
this regime, and will dominate. The LSP energy then is no longer constrained by the small 
splitting mi — m t , but by the much larger splitting m 4 ~ — m^/ — mj. If we increase the stop 
mass, the phase space for the 2-body decay begins to open up, and the distribution of LSP 
energies starts to become localized near mi — m t . But the details of this transition depend 
sensitively on the structure of the decay amplitudes. For Xi LSP, the squared amplitudes 
are directly proportional to the Xi energy. This acts to suppress the 2-body decay relative 
to the 3-body, since the latter is still capable of producing Xi with much more energy. In 
fact, the 3-body decay can contribute >10% of the branching fraction up to O(10 GeV) 
above the threshold, at least for a massless or very light x?- For the G LSP, which is 
always practically massless, the decay amplitude carries an additional momentum factor 
due to the derivative coupling of a goldstino. As a consequence, the suppression of soft 
LSP production is much stronger, and the transition region is much broader. 

For somewhat more massive neutralinos, this crossover from 3-body to 2-body decays 
occurs more rapidly as a function of m t ~. Partially this is because the spatial momentum 
available to the 2-body decay starts to open up much faster than in the nearly-massless case 
(the scaling switches from p ~ (mi — m t — m^o) to p ~ ^Jm^o(mi — m t — m~%o)), but also 
because the energy factor that appears in the matrix elements becomes bounded from 
below by m^o. However, this narrower transition to 2-body decays, where the neutralinos 
take up a smaller fraction of the available energy, is partially offset by the fact that a small 
kinetic energy carried by a massive particle can still correspond to a sizeable momentum (a 
point strongly emphasized in 21] ). Even right at the threshold, the approximately 1.3 GeV 



3 Often, it is assumed that the loop- and flavor-suppressed 2-body decay i — > ex? takes over in this 
region in the case of neutralino LSP. But the relative branching fraction between this mode and the 
flavor-conserving 4-body decay is actually highly model-dependent, and it is not difficult to make the 
flavor- violating mode subdominant [it]]. It is also possible for the stop to undergo a 3-body or 4-body 
decay through an off-shell chargino or an off-shell sfermion if cither of these is not much heavier than 
the stop. We will not explicitly consider any of these contributions in this work, though we expect that 
the off-shell chargino/sfcrmion contributions would not significantly alter our conclusions if present. 
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intrinsic width of the top quark can lead to a nontrivial momentum dispersion, for example 
at the O(10 GeV) level if the neutralino mass is 0(100 GeV). The resulting contributions 
to the missing momentum can therefore remain observably large even for very compressed 
spectra. 

The fact that there is really no way to completely squeeze out the LSP's momentum, 
even as we pass through the 2-body/3-body threshold, is advantageous for our ra^-based 
search as well as for all other searches that target the LSP's effects on the decay kinematics. 
Far above the threshold mi = m t + mLSp, the two LSP's in a stop pair production event 
carry a large amount of invisible momentum and allow m>T2 to sometimes exceed my/. At 
and below this threshold, the decays are 3-body with an on-shell W and a large fraction 
of the energy going to the LSP, again allowing for my 2 > my/. Somewhere slightly above 
the threshold, the $ T reaches a minimum, and extreme values of become much less 
common. This region is the most challenging one, but we will show that it is still often 
possible to obtain good sensitivity there. 

If we consider even lower stop masses, approaching the 4-body threshold, the rate of off- 
shell Ws rapidly climbs. Despite the added $ T , the rriT2 distribution inevitably becomes 
softer and the endpoint eventually falls below my/. This presents an absolute end to the 
utility of the search (as well as the m^ search in the Z+jets channel). Nonetheless, 
other search strategies are already being employed that cover part of this region, and we 
will describe these in subsequent sections. 

Before proceeding, we also comment on the stop decay lifetime. The decays are often 
rapid enough to be considered prompt, but not always. The simplest counterexample is 
gauge mediation with a moderately high Fsusy- For example, taking v^susy 21 100 TeV, a 
200 GeV stop can propagate more than O(mm) distance before decaying [23|, |24j. For stops 
lighter than 130 GeV, decay lengths are at least on the order of a mm unless V^Fsusy ^ 
10 TeV. Stop decays to a neutralino LSP would usually proceed much faster, unless Xi is 
dominantly singlino, in which case the lifetime can be made arbitrarily long by adjusting 
the singlino-gaugino-higgsino mixing. In this paper, we will only consider prompt decays 
and focus on traditional search strategies. This assumption can be satisfied in the bulk of 
the parameter space, except for the region m~ t < 100 GeV for a gravitino LSP, which in 
any case should have been highly visible at LEP. 

The region mi + ttilsp < m t involves an additional constraint from top quarks decaying 
into t(x\/G). However, there is a broad range of values for the stop-top-LSP couplings 
such that the branching fraction of this top decay mode is very small, while the stop decay 
remains prompt. This happens automatically in gauge mediation, and a pure bino LSP 
would be marginally safe. A pure higgsino-like LSP would probably be ruled out, but even 
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a modest mixing with a singlino state is sufficient to avoid the constraint in this case as 
well. 

B. Spin effects 

We have seen above that the momentum spectrum of the LSP in stop decay can depend 
on whether it is a neutralino or a gravitino, and (in the case of Xi) on hs mass. When we 
construct quantities like m^ 2 (^ + , 1~,$t) or m T(^-^r) i R a complete event, we must also 
pay attention to the kinematics of the leptons and neutrinos, and their correlations with 
the LSP's. This is especially true since we are looking for events that lie on the tails of 
distributions, and, as we now explain, these tails are heavily affected by the spin of the top 
along the stop decay axis. 4 

Let us first consider the case of a gravitino LSP. The couplings of the gravitino preserve 
chirality, so, for example, right-handed chirality stops couple to right-handed chirality 
top quarks. We are considering decays far below the limit where the tops are produced 
relativistically, so seemingly the top chirality and helicity are not simply related. However, 
the gravitino is relativistic, and the chirality structure of the interaction vertex completely 
governs its helicity. Right-handed chirality stops therefore always produce right-handed 
helicity gravitinos. Since the stop is spin- zero, the top quark is also forced to have right- 
handed helicity, regardless of its velocity. In other words, the top quarks are 100% polarized 
along the decay axis in right-handed chiral stop decay. In fact, this occurs even if the 
top is off-shell. For left-handed chiral stops, the same argument goes through, up to a 
subdominant contribution from the tijaiWG 4-point supercurrent coupling. 

The fact that the tops produced in chiral stop decays are polarized leads to a significant 
effect on the final state kinematics. The angular distribution of the lepton produced in 
a top quark decay is maximally correlated with the top's spin. 5 Therefore, in the parent 
stop's rest frame, the lepton tends to follow the direction of the top in right-handed decays, 
and tends to follow the direction of the gravitino in left-handed decays. There is also a 
similar bias for the neutrino, but it is less pronounced and goes in the opposite direction. 
The net effect in right-handed stop decays is that the missing energy contributions from 
the neutrino and gravitino tend to reinforce each other, and the entire invisible subsystem 



4 For a more detailed discussion of spin effects in stop decay, as well as their possible measurement at the 
LHC, see 

5 dT/<Kl(l) oc 1 + (S(t)) ■ where Cl(l) is the lepton's unit vector and (S(t)) is the top's normalized 
spin vector, both in the top's rest frame. 
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tends to be back-to-back with the lepton. This is an optimal situation for enhancing rriT2 
in the dileptonic final state and tut in the Z+jets final state in stop pair events. In contrast, 
the opposite tendency for left-handed stops tends to suppress m^i and my. For our own 
m,T2 search, we therefore expect that searches for right-handed stops decaying to gravitinos 
will be more sensitive than for left-handed stops. 

Similar effects occur with a neutralino LSP. However, there we have more choices, and the 
behavior is somewhat more complicated. First, neutralino couplings to stops come in two 
varieties: chirality-preserving gaugino couplings and chirality-flipping higgsino couplings. 
Restricting ourselves momentarily to m^o = 0, we expect exactly the same spin effects as in 
the gravitino LSP case, but with a reversed sense of chirality for a neutralino that couples 
dominantly like a higgsino. 6 There is therefore a correspondence between switching the 
LSP identity (x? ~ B) -H- (%° ~ h°) and switching the stop chirality tn <r> tz- When the 
neutralino is massive, but still couples dominantly like a gaugino or higgsino, the strength 
of the polarization depends on its velocity in the decay: fast-moving neutralinos are again 
approximately chiral and force the top to be polarized, while slow-moving neutralinos and 
their associated tops are unpolarized. The latter occurs exactly in the limit of compressed 
spectra, though the depolarization is always tempered by the presence of a nontrivial 3- 
body contribution. 

Given these observations, we can identify four distinct limiting cases for the kinematics 
of stop decay. These are: 

• i R ->tG 

• i L ->tG 

• i R -+tB or t L -> th° 

• t L ->tB or t R -> th° 

The different rriT2 distributions arising from the last two cases, for the specific example of 
220 GeV stops decaying to massless neutralinos, is illustrated in Fig. [TJ 

To fully understand the possible impact of the spin and momentum structure of the 
stop decay amplitudes, we will study the LHC sensitivity of each limiting case indepen- 
dently. While we expect that this captures most of the relevant phenomenology, a given 
spectrum can of course entail more complicated patterns of neutralino and stop mixings, 



6 There is a slight subtlety if the neutralino couples like a very pure wino. Then even a mostly-i^ will be 
forced to decay via its component. Though the coupling then technically preserves chirality, it will 
appear to flip it. In any case, our very general treatment covers even this extreme situation. 
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and interference effects could become quite nontrivial. There are nonetheless some further 
simplifications that can occur. For example, if the LSP is massless and pure gaugino, hig- 
gsino, or gravitino, but the stop is arbitrarily mixed, the contributions from the tn and ti 
components do not interfere. Also, if the LSP is couples like a mostly-bino, it preferentially 
selects out the tR component due to the larger hypercharge. 



IV. SEARCH STRATEGY 

In this section we will describe stop search strategies based on the use of dileptonic m^. 
We begin by describing the ATLAS search on which we base our own, and reproducing the 
kinematic distributions therein. This analysis was optimized for spectra far from the com- 
pressed regime and was therefore sensitive to heavier stops. Nevertheless, understanding it 
in detail gives us a chance to validate our Monte Carlo event samples, including important 
effects such as detector energy resolution. We then proceed to introduce our own version 
of this search, which we expect to yield better sensitivity in the compressed regime. 



A. The ATLAS dileptonic search 



The recent ATLAS dileptonic rriT2 search 16j is already placing limits on direct produc 



tion of right-handed stop pairs decaying as in — >• tB. The search uses the following event 
selection criteria in order to isolate the stop signal in the dileptonic channel: 

• Exactly two opposite-sign isolated leptons (I = e or fi) satisfying pr{e) > 25 GeV, 
prill) > 20 GeV, \r){l)\ < 2.5, and m{l + > 20 GeV (regardless of flavors). 

• At least two jets with pr(J) > 25 GeV and \r](j)\ < 2.5, with the leading jet satisfying 
p T (ji) > 50 GeV. 

• In same-flavor (SF) events, a Z-veto (m(7 + ,/~) ^ [71,111] GeV) and at least one 
6-tagged jet. 

• For the final signal region selection, events with rriT2 > 120 GeV. 

In order to reproduce the kinematic distributions presented in the ATLAS analysis, as 
well as for further use in our own analysis, we generate Monte Carlo samples of the signal 
and the backgrounds for SM ti, Z + /~+jets, diboson, tiW/Z and single top production. We 
do not attempt to model fake lepton backgrounds, which are unlikely to be important for 
us. The details of the Monte Carlo tools that we use, as well as our approximations of 
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^ (GeV) m T2 (GeV) 



FIG. 2: Ifjj, (left) and (right) distributions after the selection cuts of the ATLAS analysis 
have been applied, with 5 fb _1 at the 7 TeV LHC. We show ATLAS's tt prediction (red curve) and 
our own background predictions: tt (green), single-top and tiW/Z (blue), electroweak diboson 
(purple), Z + Z~+jets (orange), and total background (black solid). For we also show three 
signal points, corresponding to a massless bino-like neutralino LSP and right-handed stops with 
masses of 140 GeV (dot-dashed curve), 183 GeV (dotted curve), and 220 GeV (dashed curve). 
(Error bars are Monte Carlo statistics.) 

detector effects, are described in detail in appendix |A] For the following, we will combine 
the same-flavor and different-flavor channels into a single analysis, unlike the ATLAS search 
which keeps them separate. 

Applying the selection cuts above (except for m^), we display our distributions for 
JJjrp and rriT2 in Fig. |2j We include the tt distribution estimated by the ATLAS analysis 
for comparison. In the latter plot, we also show a small set of signal distributions, 
corresponding to t# — > tB with m^o = and mj = 140, 183, and 220 GeV. 7 The agreement 
with the ATLAS tt distributions is generally very good, with the largest statistically- 
significant discrepancies (0(20%)) occurring for the lowest two ]fi T bins. These do not 
contribute to the signal region at large rriT2- The largest resolvable discrepancy in modeling 
of the other backgrounds, which is not a very significant one, is in the / + /~+jets background. 
For m>T2 > 120 GeV we predict roughly 0.5 event, whereas ATLAS predicts 1.2 ± 0.5. In 
fact, we will soon show the / + /~+jets background to be highly subdominant for our own 

7 Though not shown, we also reproduce the shape and normalization of their example signal point, 
(m t -,m^o) = (300,50) GeV. 
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version of the m^ 2 search. 

The expected turnover near mw in the ti background is indeed present, with the endpoint 
smeared out by finite resolution effects and, to a lesser extent, by the Ws natural width. 
The falloff is nonetheless rapid, with a roughly two orders of magnitude drop between 
60 GeV and 100 GeV, and two more orders of magnitude up to 130 GeV. Until the turnover 
point, ti is the leading background, but other backgrounds that do not have as sharp a falloff 
begin to give an order-one relative contribution at very high m^- The biggest amongst 
these are tiW/Z and / + Z _ +jets. The latter can achieve large values of m^ 2 when O(10's) 
of GeV lfi T mismeasurement, often supplemented by neutrinos in 6-hadron decay, points 
opposite to a highly-boosted / + /~ system that is off of the Z resonance. 

One can see based on Fig. [2] why the ATLAS search is not very sensitive to light stops. 
The signal in the 220 GeV stop mass case exhibits a healthy S/B at high txiti-, but the 
expected event count is only 0(2) with 2011 LHC data in the region my 2 > 120 GeV. 
This feature should nonetheless become observable with increased statistics. The 140 GeV 
stop is much more difficult to observe with an cut alone, although it does achieve 
S/B ~ 1/2 for a cut near 100 GeV. The most difficult point, 183 GeV, never exceeds 
S/B ~ 1/4. 

B. Proposed analysis for compressed spectra 

Having validated our analysis setup by reproducing the ATLAS kinematic distributions, 
we attempt to optimize the dileptonic analysis for the case of light stops with masses 
in the compressed regime and below. The most obvious consequence of working in this 
region of parameter space is a reduction in the available Jfi T in signal events, and hence a 
softening of the m^ 2 distribution. As we saw above, searches for lighter stops can benefit 
greatly from lowering the final cut, and throughout the rest of our study we will utilize 
a final cut of 95 GeV. This particular cut is not always absolutely optimal for all cases, 
but it is nonetheless quite effective over a broad range of models. In particular, we have 
trained our analysis on compressed £r — > tB assuming 25 fb _1 of 8 TeV data, and expect 
this cut to be close to ideal for that case. 

In addition to reduced $ T and m,T2, we must also contend with the softer fe-jets arising in 
off-shell decays. The off-shell decays are a significant contribution for compressed spectra, 
and can yield some of the highest-mr2 events. Given this fact, it is not obviously beneficial 
to demand two hard jets in the event selection, especially since one the most efficient 
ways for the SM to produce hard dileptons with rriT2 ^ f^w at the LHC is already ti. 
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FIG. 3: rriT2 (left) and $ T /m e fi (right) distributions after our baseline selection cuts, with 5 fb _1 
at the 7 TeV LHC. We show our background predictions: ti (green), single-top and tiW/Z (blue), 
electroweak diboson (purple), Z + /~+jets (orange), and total background (black solid). We also 
show three signal points, corresponding to a massless bino-like neutralino LSP and right-handed 
stops with masses of 140 GeV (dot-dashed curve), 183 GeV (dotted curve), and 220 GeV (dashed 
curve). (Error bars are Monte Carlo statistics.) 

Consequently, we largely eliminate any requirements on jets accompanying the two leptons. 8 
A necessary exception, which helps to combat the ZZ* background, is that we continue to 
demand at least one 6-tagged jet in the SF subsample. Our baseline selection, before any 
rriT2 cut, is therefore almost the same as that of ATLAS, but with the "At least two jets..." 
requirement dropped. 

Having relaxed the jet selection cuts of the ATLAS analysis, we already gain some im- 
provement in S/B at high rriT2, as seen in Fig. [31 However, it is possible to make further 
gains. A commonly-used variable for SUSY searches, though not applied by ATLAS for 
dileptonic stops, is the ratio $ T /m e g (or $ t /St in CMS terminology). Here m e g (St) 
represents the scalar-summed px of all reconstructed objects in the event, including jets, 
leptons, and $ T . 9 We also show in Fig. [3] the distribution of $ T /m e s with our baseline se- 



8 We might also consider using lower pt cuts on the leptons themselves, and this is in fact possible thanks 
to the availability of relatively low-threshold dilepton triggers. However, by demanding high-TOT2 we are 
effectively demanding hard leptons anyway, and we do not expect a significant benefit from accepting 
softer leptons into the analysis. 

9 Other common variations are $ T /y/m c s or $ t /^JHt, the latter using the scalar pr sum over jets only. 
We find slightly better results with $ T /m e ff rather than ]£ T / -y/rn ff, but the choice of exactly how to 
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FIG. 4: $ T /m e ff distributions after application of our final m^2 > 95 GeV cut (left) and rriT2 
distributions after application of our final $ T /m e s > 1/3 cut (right), with 5 fb _1 at the 7 TeV 
LHC. We show our background predictions: ti (green), single-top and tiW/Z (blue), electroweak 
diboson (purple), Z + Z~+jets (orange), and total background (black solid). We also show three 
signal points, corresponding to a massless bino-like neutralino LSP and right-handed stops with 
masses of 140 GeV (dot-dashed curve), 183 GeV (dotted curve), and 220 GeV (dashed curve). 
(Error bars are Monte Carlo statistics.) 

lection. Some signal and background separation is evident there, especially for the 140 GeV 
stop. But in Fig. [4j we can see what happens after a 95 GeV cut. Then, even the 
difficult 183 GeV stop can be seen to achieve S/B ~ 1. This means that stop events tend 
to achieve large m^ 2 values with softer leptons and harder ]fi T relative to backgrounds such 
as ti. At the same time, we see that non-top backgrounds almost completely fall away, as 
these usually rely on smaller $ T in association with very hard dilepton systems. 

In order to capitalize on the good separation in $ T /m c s at high-m?^, we define our 
final signal region cut as $ T /m c ft > 1/3 and > 95 GeV. We can get a complementary 
picture by looking at the m T2 distribution after application of this E T /m cS cut, as is also 
shown in Fig. HJ 

Our complete event selection is then summarized as follows: 

• Exactly two opposite-sign isolated leptons (I — e or /i) satisfying pr(e) > 25 GeV, 



represent this variable is not crucial. However, $ T / V Ht is potentially dangerous because our very loose 
jet requirements can lead to a vanishing denominator. 
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Pr(fi) > 20 GeV, \rj(l)\ < 2.5, and m(l + ,l~) > 20 GeV (regardless of flavors). 

• In same-flavor (SF) events, a Z-veto (m(/ + ,/~) ^ [71,111] GeV) and at least one 
6-tagged jet. 

• For the final signal region selection, events with ^ T / m eff > 1/3 and > 95 GeV. 

For 5 fb _1 at LHC7, the dominant backgrounds are then tt (4.0 events), tW (0.6 event), 
ttW/Z (0.3 event), and ZZ* (0.2 event). For 25 fbr 1 at LHC8, the backgrounds all increase 
by approximately an order of magnitude. 

Understanding the systematic errors on these background predictions would require 
much more sophisticated detector modeling than what we have available. To proceed, we 
simply make a reasonable guess, which at a minimum will give some indication of the level 
at which these errors might affect the search sensitivity. While the nominal ATLAS analysis 
estimates an 0(1) uncertainty for the tt background at > 120 GeV, we expect that 
the errors on our own version will be smaller. In the original conference note and on the 



analysis website 32j, ATLAS reports the error on the background prediction in the range 
rriT2 = [100, 120] GeV as 15% for same-flavor and 40% for different-flavor. For our own 
95 GeV cut, we expect even smaller uncertainty, since ATLAS shows a general trend of 
decreasing error for decreasing m,T2 (as can be seen in Fig. 2 of [16J). We also expect the 
errors to decrease with larger background Monte Carlo samples passed through the full 
detector simulation, as simulation statistics are clearly still an important limiting factor 
for the ATLAS analysis. We therefore take a somewhat optimistic stance, and for our own 
analysis pick the smaller of the two errors reported for the [100, 120] GeV bin, namely 15%. 
In the next section, we will also present some of our results assuming the more pessimistic 
40% choice. 

Of course, our analysis uses different selection cuts before cutting on m^, essentially 
swapping a demand of two hard jets for a demand of somewhat large $ T /m e g. We have 
no way of assessing whether this significantly changes the error accounting. However, the 
$ T j m e s cut is not very harsh and is defined in a manner that is quite inclusive over event 
kinematics. We therefore do not expect this modification to pose a significant barrier to 
achieving a reasonable level of uncertainty. 

One place where we do have a certain degree of control over uncertainties is the theoreti- 
cal error on the background prediction. We have checked that independent 0(1) variations 
on the renormalization and factorization scales of our tt simulations, generated with MC@NL0, 
lead to variations in our predictions which are at the level of our Monte Carlo statistics. 
This is roughly 10% for events passing our final cuts. 
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V. RESULTS 



Having outlined the details of our search strategy optimized for light stops in the com- 
pressed regime, we now estimate the sensitivities that can be obtained with this analysis. 
To make contact with the existing 2011 ATLAS study, which placed exclusions for tn — > tB, 
we first describe how those exclusions might be expanded. We also show how they extend to 
the case of t# — > tG, to which there is even better sensitivity. We then show the prospects 
for the full 2012 data sample at 8 TeV. We perform our analysis for all four limiting models 
described in section [III Bl For the neutralino LSP cases, we also show a reinterpretation of 
ATLAS's t — > bxt search, which capitalizes on low-pr dileptons and further expands the 
possible coverage in the (m f ~, m^o) mass plane for very light stops. Finally we turn to the 
question of how the size of systematic errors may impact our t^ — > tB results, and we make 
projections for what may be achievable in the future. We consider both a cut-and-count 
style analysis after a high statistics run at 14 TeV, and the possibility of a precision shape 
measurement of the m>T2 spectrum. 

Below, we base our main statistical results on counting experiments using the CL$ 
procedure. As discussed above, we define our signal region (SR) using a set of baseline 
selection cuts as well as a final selection with $ T /m c s > 1/3 and > 95 GeV. Unless 
stated otherwise, we assume a 15% systematic error on the background prediction in this 
region. 

Analyses of this type also usually normalize backgrounds using control regions (CR). 
This would not be a necessary step for a simple estimate such as ours, but we include it to 
avoid misleading results. For very compressed spectra or unfavorable chirality choices, the 
$ T /m C H and rriT2 distributions can look very top-like, and a realistic analysis can poten- 
tially normalize away the signal. The actual choice of a reasonable CR would ultimately 
be determined by the experimentalists. We simply take as our CR sample the entire set 
of events satisfying our baseline selection and inverting our final selection: ^ T /m c g < 1/3 
or rriT2 < 95 GeV. (We neglect the non-top backgrounds such as diboson and / + /~+jets in 
this construction, as these can be separately normalized using alternative control regions, 
and are highly subdominant in the SR.) Effectively, our nominal statistical procedure cor- 
responds to a shape analysis with two bins. We provide complete details in appendix [Bl 

A. Prospects for 2011 data 

In Fig. [SJ we show the CL$ exclusion level estimated for the 2011 LHC data set of 5 fb _1 
at 7 TeV, for the two scenarios Ir — > tB (equivalent to ti — > th°) with massless LSP and 
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FIG. 5: Estimated exclusion level for Ir — > tB/ti, — > thP with massless LSP (left) and Ir — > tG 
(right), assuming a 2011-like data sample of 5 fb _1 at 7 TeV. The black curve and green band 
show the median and ±1<t quantile exclusions, with 15% systematic errors on the background. 
The dashed horizontal line indicates 95% CL$ exclusion. The light gray shaded region is where 
our median exclusion is better than 95%. The dark grey shaded region and vertical black lines 
are ATLAS 's observed and expected exclusions from their dileptonic rriT2 search. 



tn tG. The CLs exclusion is translated into an effective number of one-sided cr's, with 
95% exclusion corresponding to 1.64a. In addition to the median exclusion expected with a 
background-only sample, we show the the usual "green band" of la statistical uncertainty 
on the exclusion prediction. The plots also show the ATLAS 2011 exclusion using dileptonic 
my 2 , which we have re-evaluated for the gravitino case and for spectra with off-shell decays. 

For the Ir —> tB case, we can see that the potential exclusion range is greatly expanded 
on the low end, pushing into the compressed region and ruling out a large swath of the 
region with < m t . A gap nonetheless remains between 175 GeV and 210 GeV. We will 
determine the fate of this gap in the next subsection. 

For tit ^ tG, the prospects are even better, owing to the typically larger energies carried 
by the gravitino. Our median prediction is that the entire mass range mi = [105, 380] GeV 
can already be ruled out. However, the weakest point, at = 200 GeV, is just barely 
closed off by the median exclusion. It is possible that statistical fluctuations would lead to 
a weaker observed exclusion, leaving a small gap in the mass coverage. 
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FIG. 6: Estimated exclusion level for tn — > tG (left) and ii — > tG (right), assuming a 2012-like 
data sample of 25 fb^ 1 at 8 TeV. The black curve and green band show the median and ±1<7 
quantile exclusions, with 15% systematic errors on the background. The dashed horizontal black 
line indicates 95% CLg exclusion, and the dashed horizontal orange line indicates the approximate 
equivalent of "5cr" discovery level. The light gray shaded region is where our median exclusion 
is better than 95%. On the left plot, the dark gray shaded region indicates the complete range 
of ATLAS exclusions: dileptonic, /+jets, and all-hadronic. (The left edge is controlled by Z+jets, 
which we have not re- interpreted for a gravitino LSP. We expect the true exclusion to be stronger.) 
On the right plot, the dark gray shaded region indicates the ATLAS all-hadronic exclusion, which 
is likely their only search unaffected by the top quark's spin. 

B. Prospects for 2012 data 

The larger luminosity and energy of the 2012 LHC run will greatly increase the sensitivity 
of the dileptonic search. For the results presented below, we assume 25 fb _1 of data 
collected at 8 TeV. 

We begin with the gravitino LSP. In Fig. El we show how the coverage will evolve 
for tn, and now also include results for t^. The former search becomes capable of cleanly 
excluding stops between 100 GeV and 490 GeV, with no gaps. Indeed, most of the range not 
already excluded by ATLAS stop searches would exceed discovery- level significance, though 
some of this region might be independently excluded by non-dedicated SUSY searches 
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251 ] . In the ti — > tG search, we clearly see the degrading effects of left-handed top quark 



polarization. The coverage is much weaker over the entire range (note the change in vertical 
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FIG. 7: Estimated 95% exclusion region and 5a discovery region for t R ->■ tB/t L -)• th°, assuming 
a 2012-like data sample of 25 fb" 1 at 8 TeV. Our median exclusion boundary is represented by 
the solid black line with hash marks, and the ztltr quantile boundaries define the green band. 
Discoverable regions are shaded orange. We also include various existing experimental constraints. 
Low-mass LEP exclusions and the complete set of high-mass ATLAS exclusions are shaded light 
gray, with the ATLAS dileptonic rriT2 region bordered by the dotted black line. The exclusion 
boundary from CMS all-hadronic searches (inclusive razor, 6-tagged razor, and oct) is indicated 
by the dashed black line. The dotted black line with hashes shows the exclusion possible from 
the ATLAS low-p^ dilepton search for t — > bxf . Red lines indicate the boundaries between the 
different iV-body kinematic regions. 

scale). We nonetheless predict exclusion-level sensitivity in the ranges = [110, 190] GeV 
and [235, 380] GeV, and a region with discovery- level sensitivity centered at 160 GeV. 
While a gap in coverage remains, we expect this to close off as even more data comes in 
at the upgraded LHC. 10 We also emphasize that our own analysis has not been separately 

10 More generally, mixed stops may or may not exhibit this gap. For the weakest mass point, mj ~ 200 GeV, 
scanning over stop mixing angles, we estimate that the borderline case is a mostly-^ with \9f \ ~ 40°. 
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FIG. 8: A zoom-in of the compressed region of Fig. [7J using shifted coordinates. The dotted lines 
show the modification to the median exclusion boundary under ±10% variations of the signal 
cross section. 

optimized for this case, and a more careful event selection might yield better results. In 
particular, II — > tG could benefit from a harder m T2 cut. 

For Ir — > tB and ti — > th°, we extend our results into the full mass plane, shown 
in Fig. Compared to the original 2011 ATLAS dileptonic search, the possible range 
of coverage has expanded dramatically, encapsulating a large portion of both the 2-body 
and 3-body decay regions. We achieve discovery-level significance over broad parts of the 
plane, including an area between the current exclusion and the compressed region at the 
2-body/3-body interface. Our exclusion-level sensitivity extends deep into this region, just 
barely closing the gap for a massless x? an d squeezing it down to O(10) GeV width for 
higher masses. We show a more detailed zoom-in on the compressed region in Fig. [SJ There, 
we also show the effect of ±10% variations in the signal cross section, which give a rough 
indication of theoretical errors. We see that the size of the gap is not very sensitive to these 
cross section variations. As we will see in the next subsection, its extent is however highly 
dependent on systematic errors. 

We also indicate in Fig. [7J the exclusions from the full complement of ATLAS 
searches (dileptonic 16], /±jets [l5j], and all-hadronic [l4|), the various CMS all-hadronic 
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searches 17H19||. and LEP limits from a direct 3-body/4-body stop search by ALEPH 11 33] 
and from Z width constraints 34J . We supplement these limits with a re-interpretation of 
ATLAS's 2011 search for low-mass stops decaying as t — > bxt — > b(W + Xi), which capital- 
izes on low-pr dileptons 35]. This makes a substantive impact in the fully 4-body region, 
where our own strategy no longer works because the W bosons are off-shell. We have 
also studied the dileptonic portion of ATLAS's low- mass stop search based on the cluster 



transverse mass variable 36], but obtain somewhat weaker results. It would be interesting 
to make a more comprehensive study of these searches, and to better optimize them for 
direct stop decays without assuming intermediate on-shell charginos. 

By comparison, the prospects for tn — > th° and — > tB are much less optimistic. In the 
2011 data, neither our search nor ATLAS's dileptonic my 2 search are capable of placing any 
limits. The situation in 2012 is shown in Fig. |9j With rather sizeable uncertainty, we can 
exclude stops in the range [240, 380] GeV given a massless LSP. No region passes discovery- 
level significance. We note that a similar, though perhaps somewhat less severe degrading 
of sensitivity would also occur in Z+jets searches based on mr(l,^ T ), as illustrated above 
in Fig. [TJ Consequently, we do not show existing LHC exclusions for dileptonic or Z+jets 
searches for these models, but only for all-hadronic searches. However, the same LEP 
limits should continue to apply, as well as those from the ATLAS low-px dilepton and 
low cluster transverse mass searches. It would clearly be interesting here to see how the 
all-hadronic searches perform for off-shell decays, as well as how far they can be pushed 
into the compressed region using precision shape measurements (as suggested by 21]). We 
also comment that the dileptonic search proposed by 22], which focuses on distortions of 
tt spin correlation and rapidity distributions, should still apply. 



C. Effects of systematic errors and future outlook 

All of our results so far assume a 15% systematic error on the background. Over large 
portions of parameter space, the precise choice of systematic error does not significantly 
alter our estimates of exclusion or discovery potential. However, we briefly address the 
impact of these errors in the difficult compressed region for tn — > tB/ti — > th°. In Fig. [T0| 



11 In particular, we use the results of ALEPH's search for i —> blvx\ via an off-shell chargino. While we 
assume that a different diagram topology dominates, namely decay via off-shell top quark, we have found 
that the two cases have very similar decay kinematics. Chirality effects in these searches should also be 
less pronounced, since they are much more inclusive. The results presented here represent an average 
between the (very similar) constraints placed on pure a left-handed stop and a 8f = 56°, mostly-right- 
handed stop with vanishing Z coupling. 
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FIG. 9: Estimated 95% exclusion region for t R ->■ th°/t L -> tB, assuming a 2012-like data sample 
of 25 fb _1 at 8 TeV. Our median exclusion boundary is represented by the solid black line with 
hash marks, and the ±1<t quantile boundaries define the green band. We also include various 
existing experimental constraints. Low-mass LEP exclusions are shaded light gray. The ATLAS 
all-hadronic search boundary is indicated by the dotted black line. The exclusion boundary from 
CMS all-hadronic searches (inclusive razor, 6-tagged razor, and cut) is indicated by the dashed 
black line. The dotted black line with hashes shows the exclusion possible from the ATLAS low- 
Pt dilepton search for t — > bxf . Red lines indicate the boundaries between the different iV-body 
kinematic regions. 

we show how the boundaries of the exclusion change if we either optimistically assume 
vanishing systematics, or pessimistically assume 40% systematics. In the former case, the 
low-mass LSP is much more cleanly closed off, with median exclusion extending up to about 
17 GeV for any stop mass in this range. The remaining sensitivity gap is also generally 
narrower. In the 40% systematics case, the gap instead becomes much broader, already 
15 GeV wide for a massless LSP. Clearly, maintaining good control over the systematic 
errors will be crucial to narrowing the gap as much as possible, and having a good un- 
derstanding of the size of those errors will be necessary to reliably delineate the exclusion 
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FIG. 10: As in Fig. [8l a zoom-in of the compressed region of Fig. [7J using shifted coordinates. 
The systematic errors have been changed from 15% to 0% (left) and 40% (right). 

boundary in this region. 

We can also consider the improvements in sensitivity that may be achieved with a high 
statistics dataset corresponding to 100 fb" 1 at 14 TeV. To estimate this, we scale up the 
event counts used in computing the 2012 results by a factor of 15, to account for the 
enhanced luminosity as well as the increased signal and background cross sections. The 
projected sensitivity in the compressed region with these assumptions is shown in Fig. [TTJ 
It is evident from this plot that statistics were still a major limiting factor for the 2012 
analysis. The exclusion gap for a massless LSP can now be closed entirely. (At this point 
we instead have a discovery gap, of roughly 10 GeV width.) Median exclusion fully covers 
LSP's up to 28 GeV. We emphasize again that our search has not been re-optimized for 
these variations, thus higher sensitivity may well be achievable. 

Considering the power of a simple cut-and-count search using the dileptonic my 2 vari- 
able, a straightforward way to further enhance sensitivity would be to measure the detailed 
shape of the rriT2 spectrum. Our own initial studies with MC@NL0 indicate that the shape 
is quite stable to 0(1) renormalization and factorization scale variations, as might be ex- 
pected given the highly inclusive nature of our kinematic criteria. A more comprehensive 
study with higher statistics, also folding in NLO corrections to the top decay, would be in- 
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t R -> tB / t L -> th° (approx LHC14 @ 100 fb" 1 ) 




[ ' ' ' ' 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
-20 -10 10 20 30 40 
mj - m LSP - m, (GeV) 



FIG. 11: As in Fig. [8l a zoom-in of the compressed region of Fig. [7J using shifted coordinates. 
Systematic errors are again 15%, but the effective luminosity has been increased by a factor of 15, 
to roughly extrapolate to what is possible with 100 fb^ 1 at 14 TeV using our simple cut-and-count 
search. 

teresting to pursue. As a first estimate of the potential of a shape-based search, we neglect 
all possible theoretical and experimental systematic errors. We divide the distribution, 
after the $ T /m C R cut, into 5 GeV bins from 20 GeV to 125 GeV and a single overflow bin for 
events above 125 GeV. We include all backgrounds, though the top-like background (ti, tW, 
and tiW/Z) is by far the most important. In constructing the signal+background shape, 
we add the stop signal and the top-like background, which tend to have similar shapes 
at low m,T2, and rescale to the top-like background's original normalization. To build up 
probability density functions for background-only and signal+background hypotheses, we 
compute log-likelihood ratios over an ensemble of pseudo-experiments, as described in ap- 
pendix [HI A scan over the weakest (m 4 ~, m^o) points from our cut-and-count search reveals 
major improvements. For the 2012 analysis, the gap in median exclusion completely closes 
for m^o < 45 GeV. Extrapolating to 100 fb" 1 at 14 TeV, complete exclusion coverage 
extends up to m^o ~ 200 GeV. Clearly, it will be worthwhile to understand whether the 
rriT2 spectrum can be predicted accurately enough to take advantage of such an analysis. 
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In the case of a gravitino LSP, we have shown that excellent coverage can be obtained 
with the 2012 dataset using a simple cut-and-count search, if the stop is right-handed. 
But a 45 GeV gap remained for left-handed stops. Projecting to 14 TeV running with our 
nominal analysis, we expect this gap to shrink significantly, down to 5-10 GeV, and could 
likely be closed with harder cuts. The shape-based analysis for 2012 would also already 
just barely close the gap. 

Finally, £r — > th° and ti — >• tB will continue to be challenging, even at high statistics or 
with more powerful analysis procedures. For compressed and off-shell decays, the bulk of 
the signal m*r 2 shape can actually become softer than the background ti shape, though the 
high-m T2 tail can still be broader. Harder cuts at the 14 TeV LHC could potentially pick up 
the latter feature. In principle, the shape-based analysis can pick up both of these features, 
extending the 2012 exclusion for massless higgsinos down from 240 GeV to 215 GeV, and 
to masses well below 200 GeV for 14 TeV. Nonetheless, other search strategies may yield 
more immediate gains. 

VI. SUMMARY AND CONCLUSIONS 

The LHC offers our best chance to comprehensively search for the supersymmetric top 
partners that may be responsible for taming fine-tuning in the Higgs sector. Given the 
importance of this endeavor, and the current lack of signals in the most straightforward 
searches, it is becoming increasingly clear that more aggressive strategies must be em- 
ployed. In this paper, we have explored a very promising option for expanding direct stop 
searches into the poorly-covered compressed and off-shell decay regions, utilizing a simple 
modification to ATLAS's existing dileptonic rriT2 search. 

Dileptonic is perhaps the cleanest variable available for a stop search, as the high- 
rriT2 background is almost entirely dileptonic ti, and is very sharply falling beyond mw- 
Primarily, we suggest lowering the final cut used by ATLAS, which both improves 
sensitivity to lighter stops and reduces systematic errors. We also propose that this search 
be run in a more inclusive manner, essentially placing no explicit cuts on jets except for 
a 6-tag in the same-flavor channel. Non-top backgrounds can otherwise be almost entirely 
suppressed by a cut on $ T /m c fi, which also helps significantly to separate stop pairs from 
top pairs. The final separation is good enough that a simple cut-and-count analysis with 
O(10%) systematic errors is capable of making major gains in light stop coverage. 

We further propose that this search and other stop searches be applied to a broader 
class of models. The typical assumption of Ir — > tB is nominally well- motivated due to 
its simplicity, but many other possibilities exist. In particular, we have identified four 
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limiting cases with distinct exclusion or discovery prospects, encompassing six different 
model choices. These cover contributions from both stop chiralities and a broad range of 
possible LSP's, including those with gravitinos or light singlinos. 

The most promising case for the dileptonic m^ 2 search is a right-handed stop decaying 
into a gravitino LSP, which can nearly be fully excluded at low mass with no gaps using 
only the 2011 dataset, and will be definitively excluded or discovered in 2012. Left-handed 
stops are more difficult, due to a softer m,T2 spectrum arising from top polarization effects. 
Even after the 2012 run, there may be a roughly 50 GeV wide gap remaining in exclusion 
sensitivity. 

The popular tn — >■ tB case has identical phenomenology to ti — > th°. 2011 data already 
allows exclusion-level sensitivity down to 210 GeV stops with a massless LSP, with a large 
swath of off-shell decays also excluded up to just beyond m t . The 2012 data will further 
expand this range, essentially closing off light stops of any mass with a very light LSP, and 
leaving open only an O(10 GeV) wide corridor in the compressed region with more massive 
LSP's. At the same time, new regions with discovery-level sensitivity open up, including 
3-body decays with LSP's up to 150 GeV, and an interesting semi-compressed region where 
the stop undergoes a 2-body decay with only 20~40 GeV of available kinetic energy. 

An alternative set of models features tn — > th°, or equivalently t^ —> tB. These tend 
to be the most difficult to search for using dileptonic m^, as they have the least favorable 
combination of spin and momentum- scaling effects in their decays. We estimate that these 
cases are unobservable in 2011 searches, and only marginally visible in 2012 searches. 
We emphasize that similar degrading effects also afflict /+jets searches based on m T (I,^ T ), 
though they should be somewhat less pronounced. While we expect dileptonic rriT2 to play 
a more limited role in searches for these models, it is not clear which alternative search 
would offer the best sensitivity, or whether the compressed regime can be fully covered for 
any LSP mass. One of the best remaining options may be the precision measurement of ti 
spin correlations and rapidity distributions, as proposed by [22]. 

Our work can point the way to possible future improvements and extensions. An obvious 
question is whether the small gap in coverage for compressed t^ — > tB can ever be fully 
closed. Significant progress can certainly be made by the 14 TeV LHC, just exploiting 
the much higher statistics. In order to obtain higher sensitivity still, one may attempt a 
more ambitious rriT2 shape analysis rather than a cut-and-count approach. Our preliminary 
results suggest that a shape analysis is capable of closing the exclusion gap for neutralinos 
up to 45 GeV with 2012 data, and potentially up to 200 GeV with 14 TeV data. The 
other models that are difficult to cover with the simple cut-and-count search, because they 
involve left-handed top quarks, also become much more feasible. 



28 



We have seen that while our method is extremely efficient in covering the 3-body stop 
decay region, it gives no coverage when the W bosons are forced off-shell in the fully 4-body 
stop decay. We have indicated that a re-interpretation of the low-mass t — > b\i searches 
by ATLAS already leads to some progress in covering the 4-body region. It would be 
interesting to explore how far such searches can be pushed. Conversely we expect that our 
rriT2 approach can lead to progress in searching for t — > b\i in regions of parameter space 
where the W boson is on-shell. More generally, searches for various other SUSY spectra or 
non-SUSY models can benefit, in particular searches for fermionic top partners. 

Finally, we comment that all of our signal modeling has been done at leading-order, 
supplemented by parton showers and in some cases shower/matrix-element matching. To 
conclusively understand what is happening for compressed spectra, where we have seen 
that the kinematics can be subtle, it would be very useful to have access to more precise 
predictions for rriT2 distributions modeled at NLO in both production and decay, and to 
systematically understand the uncertainties on these predictions. Automated NLO cor- 
rections on the production side should soon be possible with the MadGolem program 37]. 
Kinematic distributions for first and second generation squarks, complete through NLO in 
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production and decay, are also already being studied 

Despite the unprecedented center of mass energy and integrated luminosity of the LHC 
so far, we are still at the early stages of the LHC era. While searches will continue to probe 
ever- increasing mass scales, it is important to continue working to understand whether we 
are already copiously producing relatively light new particles which are nonetheless buried 
amidst the huge Standard Model backgrounds. This may in fact be the case for supersym- 
metry. We have shown in this paper that even some of the most difficult supersymmetric 
spectra are not hopelessly buried, and we have pointed to future improvements that can 
be achieved using the huge statistics expected over the lifetime of the LHC and with more 
sophisticated analysis techniques. 



Note added: While this work was being completed, CMS announced their version of the 
direct stop search using Z+jets and m-rfl,^), utilizing almost 10 fb _1 of data at 8 TeV 



A major difference between CMS and ATLAS is that the former uses PYTHIA6 40| for its 



SUSY simulations, losing all spin effects, whereas the latter uses Herwig++ [4l|, which has 
a complete treatment of spin. Indeed, the CMS results are weaker than what might be 
expected by extrapolating from ATLAS 's results, and CMS has pointed out that this is 
because the top quarks in ATLAS 's simulations are polarized. 
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Appendix A: Simulation 



We generate our ti background samples using the MC@NL0 event generator 42|, |43j, nor- 
malizing the cross section to 165 pb (230 pb) for LHC7 (LHC8) (see, e.g., [44J and references 
therein). For all of our other samples, including additional ti samples for cross checks, we 
use MadGraph5 vl . 4 . 7 [45J interfaced with PYTHIA |40J] with default settings. For sig- 
nals with a neutralino LSP, we use the MSSM model built into MadGraph5. For signals 
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47|. In all 



with a gravitino LSP, we use the publicly available stopnlsp UFO model 
simulations, we set m t = 172.5 GeV and T t = 1.3 GeV. 

Signal samples are treated as simple 2—7-2 processes followed by 3- or 4-body de- 
cays into bl + uxi or bl + uG, with cross sections normalized to the NLO+NLL predictions 
for 7 TeV [4s| . For 8 TeV running, we extrapolate using the leading-order ratios from 
MadGraph5, which agree with NLO ratios derived from PR0SPIN0 49j. We set stop widths 
to 0.01 GeV. All other SUSY particles are decoupled. For a small subset of models, we 
have generated samples matched up to one extra jet in production, and found no signif- 
icant differences in the relevant kinematic distributions. This gives some indication that 
our claimed final signal efficiencies are under good theoretical control, at least on the pro- 
duction side. It would also be useful to have a more systematic treatment of stop decays 
beyond leading-order. We note that some degree of matrix element and shower matching 
in the decay is automatically provided by PYTHIA 40j, when the top quark appears as an 
on-shell particle in the MadGraph5 event record. 

For backgrounds other than ti, we often match up to one or two extra jets, using the 
5-flavor cone-MLM prescription 50| with pt = 20 GeV and AR = 0.3. We apply A-factors 
to rescale the cross sections to their NLO values (which can be found, for example, in 161]) . 
The background samples include: 

• tW matched up to one extra jet (A = 1.15). 

• tiW and tiZ (K = 1.3). 
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• W + W matched up to two extra jets (K = 1.4). 

• WZ and ZZ (K = 1.4). 

• / + / _ +jets matched up to two extra jets (K = 1.3). 

After showering and hadronization, we reconstruct events in the manner described 
in [3], using Fast Jet v2.4.2 for jet reconstruction 51|. and smearing final object en- 
ergies and directions by-hand without an explicit detector model. Jets are defined with 



the anti-Zc^ algorithm with radius R = 0.4 52j. We smear the jet energies by a relative 



fraction (5 GeV)/E © (0.5 GeV 1/2 )/VE © 0.04, and r//0 individually by 0.025. We smear 
electron energies by 0.02, and muon energies by (0.1) x y/E/TW. We define $ T as the 
2- vector momentum sum of all final-state invisible particles (neutrinos and LSPs). We 
smear the individual x and y components of this vector by (0.7 GeV 1 / 2 ) x^Hx, where 
is the scalar-summed px of all visible particles with \r]\ < 5. 12 We calibrated this coefficient 
within our own simulation such that the spectrum for ti matches ATLAS 's simulation, 
which in turn agrees well with the data when added to the other backgrounds. The coef- 
ficient also matches ATLAS's most recent measurement of E T resolution, based on 2011 
data [53|. (Note that the resolution degrades with pileup, but that subtraction methods 
should actually be capable of stabilizing the coefficient at 0.5 GeV 1 / 2 .) To match the nor- 
malization of the ATLAS ti simulation, we reweight all events by an efficiency factor of 0.8. 
This effectively incorporates effects such as inefficiencies in isolated lepton identification. 

We identify a jet's flavor by matching in r]-<f) space to bottom- and charm-flavored 
hadrons in the event record. The assigned flavor is then the heaviest one found. We apply 
6-tagging efficiencies of 0.6, 0.1, and 0.01 for 6-jets, c-jets, and unflavored-jets, respectively. 
We have found that our relative rates of 6-tagged and untagged ti events are in good 
agreement with ATLAS. 



Appendix B: Details of Statistical Procedures 

Our starting point is the signal and aggregated background cross section predictions in 
the control region (CR) and signal region (SR): o^cr, 0"s,sr, ctb,cr, 0"b,cr- From these, we 



12 Since we include all particles, the amount of smearing is sensitive to the shower and underlying event 
models. This is not a major issue as long as all samples use the same models. However, the MCQNLO 
Ht spectrum for it comes out about 10% harder than the MadGraph5 Ht spectrum for 2 — > 2 it. 
Nonetheless, we have checked that a 10% rescaling of Ht has only a minor impact on the final smeared 
witi distributions. 
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define relative efficiencies to fall into the SR, 



e s = — , e B = • (Bl) 

OS.CR °B,CR 

We also define the signal fraction in the CR, 

/or - (B2) 

^BjCR 

For a given "observed" CR count iV CR , we can then predict what would be expected in the 
SR for the background-only hypothesis and the signal-plus-background hypothesis: 

B SR = e B N CR , (S + B) SR = eB + e ? /cR iVcR- (B3) 

1 + JCR 

Our test statistic is the "observed" SR count iVg R . Note that this method automatically 
loses discrimination power if eg = e#, in which case a realistic experiment would indeed 
normalize-away the signal. 

We determine the exclusion potential of the signal using the CL$ method 54(. Based 
on £?sr and (S + B)sr, we run Poissonian pseudo-experiments to build up the discrete 
probability density functions (p.d.f.'s) in N§r. We fix Nqr to its expectation value for 
background-only for a given luminosity. (The effects of introducing Poisson fluctuations 
into iVcR, or of instead setting it according to the (S + B) expectation, are small and are 
not considered.) The fact that the quantiles are discrete leads to a minor ambiguity over 
which bins to include in the CL$ calculation. We always err on the conservative side, and 
pick the weaker of the two choices. 

In order to account for systematic errors on the background modeling, we also allow 
ee to fluctuate in our pseudoexperiments. We typically assume 15% gaussian error. We 
do not apply a systematic error on the signal efficiency or / CR . (eg is typically much less 
affected by instrumentation effects than es, and those effects are correlated between the 
two e's. fcR is affected by theoretical uncertainties in the total cross section, which we 
briefly discuss in section IV Bl ) 

As a possible extension of our search, we also study a full shape-based analysis with 
vanishing systematic errors. These compare binned m^2 distributions (after the $ T /m e g 
cut) of pure background and signal-plus-background. While we include all backgrounds in 
this analysis, we separately rescale the signal added to the top-like background (ti, tW, and 
tiW/Z) to match the normalization of the top-like background alone. Our test statistic is 
the log-likelihood ratio of B and S + B hypotheses over all bins (essentially the sum of the 
measured bin counts, weighted individually by (S + B)i/Bi — 1). We construct the p.d.f.'s 
in this test statistic using sets of pseudo-experiments distributed according to either B or 
S + B. 



32 



Some of our results claim very high exclusion potential, for example in excess of 5cr. 
Estimating CLs in these cases becomes impractical using a pseudo-experiment approach. 
In practice, we use 10 5 pseudo-experiments to generate each p.d.f. When our CLs com- 
putation involves integrals that contain less than 5 events, we switch to a gaussian model. 
When they contain between 5 and 50 events, we average the exclusions (expressed in ex's) 
from the two procedures. 

For setting discovery thresholds, we determine whether the median of the S + B p.d.f. is 
separated by more than 5a from the median of the B p.d.f., measured in standard deviations 
of the latter. We always use the gaussian approximation for this calculation. Usually, the 
threshold for 5a discovery has a close correspondence with a particular CL S exclusion level, 
and we have exploited this correspondence in constructing the discovery line in Fig. |6j 

Note that our treatment simply sums over same-flavor and different-flavor channels. A 
more comprehensive analysis would treat them independently. In particular, some model 
points with off-shell decays may have low efficiency for 6-jet reconstruction, which would 
deplete the same-flavor channel's S/B due to our 6-tag requirement. 
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